Nuclear halo structure and pseudo-spin symmetry 
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Nuclear halo structure and conservation of relativistic symmetry are studied within the framework 
of the relativistic Hartree-Fock-Bogoliubov (RHFB) theory. Giant halos as well as ordinary ones are 
found in Cerium isotopes close to the neutron drip line. Bridged by T = channel, the conservation 
of pseudo-spin symmetry (PSS) plays an essential role in stabilizing the neutron halo structures. The 
Fock terms, especially the p-tensor couplings, not only play significant role in the PSS conservation 
but also present substantial contributions to the T = channel, from which is well demonstrated 
the necessity of Fock terms. 
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Since the neutron halo was first proven to exist in n Li 
[l[ , unexpected exotic modes have intensively challenged 
our understanding of exotic nuclei with extreme neutron- 
to-proton ratio, that play an essential role in the evolu- 
tion of the universe. As the typical exotic mode, nuclear 
halo (see Refs. and references therein) - the ex- 

tremely diffuse matter - may strongly enhance the reac- 
tion cross section, which is of special significance for the 
element synthesis in astrophysics as well as in the dis- 
covery of new superheavies. Accompanying nuclear halo 
occurrence, shell quenching 0, 0] has been found when 
approaching these exotic regions, e.g., the N = 8 shell in 
Li. 

The nucleon-nucleon (NN) interaction is originally due 
to meson exchange processes as predicted by Yukawa Q . 
Within such approach the nuclear binding is achieved 
mainly by the equilibrium between the scalar (phe- 
nomenological a) and vector (w) meson fields Q whereas 
the tensor forces due to p and 7r meson exchanges take 
part in the shell structure evolution flQll . As a relativis- 
tic symmetry in the Dirac equation , the pseudo-spin 
symmetry (PSS) [H, EH is an important general fea- 
ture in the nuclear energy spectra. Rs origin is due to 
a Lorentz scalar potential (a field) and a Lorentz vector 
potential (time component of uj field) equal in magnitude, 
but opposite in sign [lij]. In fact, the competition of large 
scalar and vector fields in nuclei explains naturally the 
spin-orbit potential [l5| . which makes the exploration in 
exotic regions more reliable. 

In exotic nuclei the valence neutrons or protons are 
loosely bound and their coupling with the continuum 
becomes important. In terms of the Bogoliubov quasi- 
particleSjthe relativistic Hartrec-Bogoliubov (RHB) the- 
ory [Ra-Qjl provides a unified and self-consistent descrip- 



tion of both mean field and pairing correlations and auto- 
matically takes continuum effects into account. Besides 
11 Li [18|], the halo phenomena have been predicted by 
the RHB theory in Nc 
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topes. Giant halo structures in Zr isotopes have also been 
studied within the RHB framework |23| as well as in the 
non- relativistic Hartree-Fock-Bogoliubov method [24| • 

While limited by the Hartrce approach, important in- 
gredients such as the spin- dependent tensor forces are 
missing in the RHB theory. In the density dependent 
relativistic Hartrce- Fock (DDRHF) theory [H,|26j], the 
tensor forces due to tt and p meson exchanges can be 
naturally taken into account and have brought significant 
improvements on the consistent description of shell evo- 
lution [27} and appropriate conservation of PSS [26l [28|. 
In this work, we use Cerium isotopes to study the nu- 
clear halo phenomenon and relevant conservation of PSS 
within the relativistic Hartree-Fock-Bogoliubov (RHFB) 
theory (29[, an extension of DDRHF. For Cerium iso- 
topes the proton number Z = 58 is closely related to 
the pseudo-spin partner states 7rl/ (7rlg 7 / 2 and 7r2d 5 / 2 ), 
which represent well conserved pseudo-spin symmetry ei- 
ther from the experimental data l30l or from the theoret- 
ical calculations 
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As we will see, the corresponding 
PSS conservation is essential for the neutron shell effects 
when approaching the neutron drip line. 

In Fig. [1] we show nuclear matter distributions (left 
panels) and neutron canonical single-particle configura- 
tions (right panels) for Cerium isotopes close to the drip 
line. The neutron drip line calculated here with PKA1 
[H[ is N = 140 whereas the calculations with PKOl [H[ 
and DD-ME2 [Hi predict a shorter one as N = 126. 
As shown in Figfljt, the neutron densities become more 
and more diffuse after the isotopes 186 Ce (N = 128), a 
direct and distinct evidence of halo occurrence. From 
Fig. [TJd we can see that such extremely extensive mat- 
ter distribution, e.g. in 198 Ce, is mainly due to the low-^ 
states, namely the halo orbits f4si/2, ^3<i 5 / 2 and ^3^3/2- 
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FIG. 1: (color online) (a) Neutron and proton densities, (b) the relative contributions of different orbits to the full neutron 
density in 198 Ce, (c) neutron canonical single-particle energies occupation probability (in x-error bars) and Fermi energy Ef 
(in open circles), and (d) neutron numbers filling in the halo orbits 4si/2, 3d 5 /2, and 3d 3 / 2 (TViiaio), an d the ones lying beyond 
the spheres with the radii r = 10, 11, 12, 13, 14, 15, 16fm, respectively A r r >iofm, A r r>iifm, TVr>i2f m , TVr>i3fm, A r r>i4fm, TV r> i 5fm 
and iV r> i6f m . The results are calculated by RHFB with PKA1 (2(| plus the Gogny pairing force D1S [32] ■ The spherical box 
radius is adopted as J? m ax = 28fm. 



Seen from the occupations of the halo orbits TVhaio m 
FigUJi, the halos in 186 Ce, 188 Ce and 190 Ce are ordinary, 
while 192 Ce, 194 Ce, 196 Ce and 198 Ce presumably have gi- 
ant halos because more than two neutrons are occupying 
the halo orbits. Similar conclusions can be also obtained 
from the neutron numbers lying beyond the sphere with 
the radius r = lOfm (A r r> iof m in Fig. [Hi), which is large 
enough (the neutron matter radius r n = 6.2fm in 198 Ce) 
for halos. Even extending to r = 16fm, which is sometime 
taken as the radial cut-off in the calculations of stable nu- 
clei, there are still some substantial amount of neutrons 
lying beyond this sphere in the isotopes from 192 Ce to 
198 Ce. 

In FigJTJ; we find that the halo orbits (i/As 1 / 2 , v3d 5 / 2 
and f 3^3/2) are located around the particle continuum 
threshold where they are gradually occupied. For the iso- 
topes beyond 184 Ce (TV = 126), the Fermi levels (in open 
circles) approach the continuum threshold rather closely 
such that the stability of these halo isotopes becomes 
sensitive to pairing effects. Nearby the low-? states, we 
find the high-/ states v2g 9 / 2 and v2g 7 / 2 - Because of the 
relatively large centrifugal barrier for g-orbits they do 
not contribute much to the diffuse neutron distributions. 



Nonetheless, the existence of the high-Z states nearby halo 
orbits is still particulary significant, because it leads to a 
rather high level density around the Fermi surface, and 
evidently the pairing effects are enhanced to stabilize the 
halo isotopes. 

Evidence for the existence of a halo can also be found 
by studying the systematic behavior of nuclear bulk prop- 
erties such as radii. In Fig. [5] we show the isospin depen- 
dence of the neutron skin thickness (r„ — r p ) calculated in 
RHFB theory using the parameter set PKA1 for Ca, Ni, 
Zr, Sn and Ce isotopes. With respect to the behavior in 
the stable region (shown with the dashed lines), contin- 
uously growing deviations are found in the chains of Ca, 
Zr and Ce until the neutron drip line. This can be con- 
sidered as evidence of a halo. Despite the deviations in 
the mid-region, Ni as well as Sn show an identical isospin 
dependence in both stable and neutron drip line regions, 
which may indicate only a neutron skin since the growth 
of a halo is interrupted. Compared with Ni, Sn isotopes 
show a much weaker skin effect. 

Similar systcmatics as those of Fig. 2 are also found 
in RHB calculations for the chains of Ca and Zr (2^, H^) . 
For the Ce isotopes the situation becomes quite differ- 
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FIG. 2: (color online) Neutron skin thickness r n — r p (fm) for 
Ca, Zr, Ni, Sn and Ce as a function of the isospin (N — Z). 
r„ and r p are respectively the neutron and proton root mean 
square radii. The results are calculated by RHFB with PKA1 
[29 ] plus the paring force D1S [3^ |. 
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ent. RHFB calculations with the parameter set PKAl 
(see Figs. [T]and[2]) show clear evidence for the existence 
of halo structures in the Ce isotopes, while in the calcula- 
tions of RHFB with PKOl and RHB with DD-ME2 
[33j . the isotopic chain ends at N = 126, before the halo 
occurrence predicted by RHFB with PKAl. This devi- 
ation between models can be preliminarily interpreted 
by the shell structure evolution in Fig. [3^, where much 
stronger shell effects are provided with PKOl and DD- 
ME2 than with PKAl. As shown in FigHJ:, the neutron 
shell gap (N = 126) between ^Hi3/2 an d v2g$/ 2 states 
is close to the particle continuum threshold, which might 
essentially influence the stability of the drip line isotopes. 

In Fig. [3jo, the pseudo-spin orbital splitting A£ t1 j 
given by PKAl shows an isospin dependence consistent 
with the neutron shell evolution, whereas in the results of 
PKOl and DD-ME2 such consistency is destroyed with 
the violation of PSS on the pseudo-spin partner states 
7rl/. To clarify this consistency between neutron shell 
evolution and proton PSS conservation we present in Fig. 
S] the two-body interaction matrix elements V a b calcu- 
lated with PKAl and responsible for the coupling be- 
tween the proton (a: 7r2d 5 /2 (filled symbols) and Trlg7/ 2 
(open symbols)) and neutron valence orbits (b: 1/2/7/2, 
" 2/ 5 /2, ^3p 3 /2 and v3p t / 2 (FigHa), and vlh 9 / 2 ,vli 13 / 2 
and v2gQ/ 2 (FigHJa)). It is found that the neutron orbits 
with nodes (i/2f 7 / 2 , ^2/ 5/2 , ^3p 3 / 2 , v3pi/ 2 and ^2g 9 / 2 ) 
show a stronger coupling with the proton n2d 5 / 2 than 
with 7rlg 7 / 2 orbit, while those without node (vlhg/2 and 
^1*13/2) exhibit the opposite trend. 

From 142 Ce (N = 84) to 148 Ce (N = 90) the valence 
neutrons are mainly filling the v2f 7 / 2 orbit, which leads 
to a corresponding recovery of PSS, as shown in Fig. 
[3}d. After 148 Ce the orbit v\hg/ 2: which has a stronger 
coupling with the 7rl 177/2 than with the Tr2dc,/ 2 orbit, 

starts to be gradually occupied. The PSS on nlf is 
still well preserved in a fairly long range from N = 90 



FIG. 3: (a) Neutron shell gap at N = 126 (AE N=12e , = 
£v2 Ss) /, — £vii 13 /,) and (b) proton pseudo-spin orbital split- 
ting (AE nl j: — E n 2d 5/2 ~ E 7r i g7/2 ) as functions of neutron 
number N for Cerium isotopes. The results are calculated by 
RHFB with PKAl [H, PKOl [2I, and by RHB with DD- 
ME2 [33| . The Gogny force D1S |32f | is adopted in the pairing 
channel. 

to N = 110, due to the equilibrium between the states 
with nodes (^2/ 7 / 2 , u2f 5 / 2 , ^3p 3 / 2 , v3pi/ 2 ) and the one 
without node (i/lhg/ 2 ). From N = 112 to 126, valence 
neutrons are filling the ^l«i3/2 level, as the orbits below 
are nearly fully occupied. An increase of AE 1Tl j is there- 
fore found (see Fig. \3jp) . When approaching the neutron 
drip line, the state v2g^j 2 starts to be occupied, and the 
levels below that are fully occupied. It leads as a result 
to the continuous recovery of PSS after 184 Ce. In fact 
not only the ^2g 9 / 2 state, but also the other valence or- 
bits beyond N = 126 do play a positive role in the PSS 
recovery. From Fig. [3] and Fig. [4] one can find that the 
nodal structure shows substantial effects in determining 
the NN interaction strength and further play an essential 
role in recovering the PSS. In fact we also found similar 
nodal effects as shown in Fig. [4] from the calculations 
with PKOl and DD-ME2 while the unphysically large 
gap between 7rlg 7 / 2 and ir2d 5 / 2 breaks the consistency 
(see Fig. [3]) because the valence protons can only occupy 
the 7rl<77/2 orbit. 

Comparing PKAl to PKOl and DD-ME2, one finds 
that the deviations on the neutron shell effects at N = 
126 in Fig. [3^, are mainly due to the proton influence, 
i.e., the PSS conservation. In Fig. 2}d it is seen that the 
protons occupying the 71-137/2 level tend to enhance the 
neutron shell effects (AE]\[ = i 2 e), whereas much weaker 
effects are contributed by particles in the n2d 5 / 2 orbit. 

Because of the serious violation of PSS on irlf given by 
PKOl and DD-ME2, the valence protons occupy only 
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FIG. 4: (color online) The two-body interaction matrix el- 
ements V a b- The notations are a — ir2d 5 / 2 (filled symbols), 
7T137/2 (open symbols); and 6 = v3p 1/2 , v3p i/2 , v2f 5/2 , u2f 7/2 
(a), and b — vlhg/ 2 , ^1*13/2, v ^99/2 (b). The results corre- 
spond to RHFB with PKA1 plus the pairing force Gogny 
D1S. 

the 7rl(77/2 state and this greatly enlarges neutron shell 
effects. It results in the interruption of the isotopic chain 
at N = 126. In contrast the PSS is properly conserved 
by PKA1 such that the pseudo-spin partner states are si- 
multaneously occupied and the occupations are changed 
consistently. It explains well the consistency between 
neutron shell effects and proton PSS conservation in Fig. 
[3j Although some minor differences may exist on the 
neutron side between models, this would not bring any 
substantial changes. 

Since the PSS conservation is tightly related with neu- 
tron shell effects at N = 126, it is also essential for the 
stability of Cerium halo structures. Such consistency is 
bridged by the coupling between valence neutrons and 
protons, i.e., the NN interaction in T = channel. In 
Table U it is seen that Fock terms present significant con- 
tributions to the T = NN interaction, about 30% of the 
total in most cases. In the T = channel, Fock terms are 
completely managed by isovector mesons, mainly (~2/3) 
by the tensor p, which can not be efficiently taken into 
account by the Hartree approach. From Table Q] one can 
find that the tensor p plays a significant role not only in 
the conservation of PSS [26J but also in bridging neutron 



halo structure and proton PSS recovery. 

In summary, within the relativistic Hartree-Fock- 
Bogoliubov (RHFB) theory with density dependent 

TABLE I: Two-body interaction matrix elements V a b (MeV) 
between proton and neutron valence orbits in 198 Ce, as well 
as the ratios of contributions from different coupling chan- 
nels. The four last columns give detailed Fock contributions, 
respectively from p- vector (p-V), p- vector-tensor (p-VT), p- 
tensor (p-T) and 7r-pseudo-vector (7T-PV) couplings. The re- 
sults correspond to RHFB with PKA1 plus the pairing force 
Gogny D1S. 





V a b Va b Vab 


p-V p-VT p-T tt-PV 


b 


a = Tv2d 5 / 2 




-0.257 63.2% 36.8% 


9.7% -4.0% 29.1% 2.1% 


v2gg /2 
v4s 1/2 
v3d 5/2 
u3d 3/2 
v2g 7l2 


-0.217 64.6% 35.4% 
-0.097 71.8% 28.2% 
-0.158 68.4% 31.6% 
-0.156 67.7% 32.3% 
-0.159 56.8% 43.2% 


8.9% -3.5% 27.0% 3.0% 
7.1% -4.4% 21.5% 3.9% 
8.9% -5.4% 25.7% 2.5% 
7.9% -5.4% 22.4% 7.4% 
9.3% -8.1% 30.8% 11.3% 


b 




fl»13/2 


-0.386 63.5% 36.5% 


8.3% -4.4% 22.9% 9.8% 


^99/2 
vAs 1/2 
v3d 5/2 
v3d 3/2 
v2g 7/2 


-0.157 64.1% 35.9% 
-0.052 68.1% 31.9% 
-0.078 67.8% 32.2% 
-0.074 70.2% 29.8% 
-0.138 69.0% 31.0% 


8.1% -6.3% 25.3% 8.8% 
8.2% -5.0% 25.2% 3.5% 
7.9% -5.7% 24.5% 5.6% 
8.4% -5.4% 25.9% 0.9% 
9.2% -5.0% 28.0% -1.2% 



meson-nucleon couplings, we studied nuclear halo phe- 
nomena occurring in Cerium isotopes and the relevant 
PSS conservation, and the role of Fock terms therein. 
Giant halos as well as ordinary ones are found in the 
drip line Ce isotopes. We also found that the stability 
of neutron halo structures is tightly related to the PSS 
conservation on the proton side. The Fock terms, mainly 
the p-tensor couplings, present substantial contributions 
to the NN interaction in T = channel, which accounts 
for the consistency between neutron halo structures and 
proton PSS conservation. In addition, the necessity of 
Fock terms, especially the tensor p, is well demonstrated 
as such effects can not be efficiently obtained with the 
Hartree approach. 

The author W.H.L. would like to thank Prof. H. 
Sagawa for fruitful discussions. This work was supported 
by the Alexander von Humboldt Foundation, and Major 
State 973 Program 2007CB815000, as well as the Na- 
tional Natural Science Foundation of China under Grants 
No. 10435010, No. 10775004, and No. 10221003, and 
by the U.S.DOE grants DE-FG02-08ER41533 and DE- 
FC02-O7ER41588 (UNEDF, SciDAC-2), and by the Re- 
search Corporation. 



5 



[1] I. Tanihata, H. Hamagaki, O. Hashimoto, Y. Shida, 
N. Yoshikawa, K. Sugimoto, O. Yamakawa, 
T. Kobayashi, and N. Takahashi, Phys. Rev. Lett. 
55, 2676 (1985). 

[2] J. E. Howard, H. R. Dullin, and M. Horanyi, Phys. Rev. 
Lett. 84, 3244 (2000). 

[3] P. M. andl. A. Sulai, A. C. C. Villari, J. A. Alcantara- 
Nunez, R. Alves-Conde, K. Bailey, G. F. Drake, 
M. Dubois, C. Eleon, G. Gaubert, R. J. Holt, R. F. 
Janssens, N. Lecesne, Z.-T. Lu, T. P. O'Connor, M.- 
G. Saint-Laurent, J.-C. Thomas, and L.-B. Wang, Phys. 
Rev. Lett. 99, 252501 (2007). 

[4] V. Rotival and T. Duguet, Phys. Rev. 79, 054308 (2009). 

[5] V. Rotival, K. Bennaceur, and T. Duguet, Phys. Rev. 79, 
054309 (2009). 

[6] J. Dobaczewski, I. Hamamoto, W. Nazarewicz, and J. A. 

Sheikh, Phys. Rev. Lett. 72, 981 (1994). 
[7] B. Chen, J. Dobaczewski, K. L. Kratz, K. Langanke, 

B. Pfeiffer, F. K. Thielemann, and P. Vogel, Phys. Lett. 

B 355, 37 (1995). 
[8] H. Yukawa, Proc. Phys. Math. Soc. Japan 17, 48 (1935). 
[9] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 

(1986). 

[10] T. Otsuka, T. Suzuki, R. Fujimoto, H. Grawe, and 
Y. Akaishi, Phys. Rev. Lett. 95, 232502 (2005). 

[11] J. N. Ginocchio, Phys. Rep. 414, 165 (2005). 

[12] A. Arima, M. Harvey, and K. Shimizu, Phys. Lett. B 30, 
517 (1969). 

[13] K. Hecht and A. Adler, Nucl. Phys. A 137, 129 (1969). 
[14] J. N. Ginocchio, Phys. Rev. Lett. 78, 436 (1997). 
[15] J. Meng, H. Toki, S. G. Zhou, S. Q. Zhang, W. H. Long, 
and L. S. Ceng, Prog. Part. Nucl. Phys. 57, 470 (2006). 
[16] J. Meng, Nucl. Phys. A 635, 3 (1998). 



[17] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and 

P. Ring, Phys. Rep. 409, 101 (2005). 
[18] J. Meng and P. Ring, Phys. Rev. Lett. 77, 3963 (1996). 
[19] W. Poschl, D. Vretenar, G. A. Lalazissis, and P. Ring, 

Phys. Rev. Lett. 79, 3841 (1997). 
[20] G. A. Lalazissis, D. Vretenar, W. Poschl, and P. Ring, 

Nucl. Phys. A 632, 363 (1998). 
[21] J. Meng, I. Tanihata, and S. Yamaji, Phys. Lett. B 419, 

1 (1998). 

[22] J. Meng, H. Toki, J. Y. Zeng, S. Q. Zhang, and S. G. 

Zhou, Phys. Rev. C 65, 041302 (2002). 
[23] J. Meng and P. Ring, Phys. Rev. Lett. 80, 460 (1998). 
[24] M. Grasso, S. Yoshida, N. Sandulescu, and N. V. Giai, 

Phys. Rev. C 74, 064317 (2006). 
[25] W. H. Long, N. V. Giai, and J. Meng, Phys. Lett. B 640, 

150 (2006). 

[26] W. H. Long, H. Sagawa, N. V. Giai, and J. Meng, Phys. 

Rev. C 76, 034314 (2007). 
[27] W. H. Long, H. Sagawa, J. Meng, and N. V. Giai, Euro- 
physics Letters 82, 12001 (2008). 
[28] W. H. Long, H. Sagawa, J. Meng, and N. V. Giai, Phys. 

Lett. B 639, 242 (2006). 
[29] W. H. Long, P. R ing, N. V. Giai, and J. Meng(2009), 

|nucl-th/0812.1103| 
[30] Y. Nagai, J. Styczen, M. Piiparinen, and P. Kleinheinz, 

Phys. Rev. Lett. 47, 1259 (1981). 
[31] W. H. Long, T. Nakatsukasa, H. Sagawa, J. Meng, 

H. Nakada, and Y. Zhang, Phys. Lett. B 680, 428 (2009). 
[32] J. F. Berger, M. Girod, and D. Gogny, Nucl. Phys. A 

428, 23 (1984). 
[33] G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, 

Phys. Rev. C 71, 024312 (2005). 



